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CnI ■ Abstract 

We consider models of bosons on curved 3+1 dimensional space-time embedded 
in a higher dimensional flat ambient space. We propose to derive (rather than 
postulate) equations of motions by assuming that a standard Klein-Gordon field on 
ambient space is restricted to space-time by a strong confining potential. This leads 
^ ^ to a modified Klein-Gordon equation on space-time which includes, in addition 

to the standard terms, a term with a so-called induced potential which depends 
5^ , on intrinsic- and extrinsic curvature of the embedded space-time but not on the 

details of the confining potential. We compute this induced potential for natural, 
simple embeddings of Schwarzschild- and Robertson- Walker space-times. We also 
. discuss possible observable implications of our results and, in particular, propose and 

I study an extension of a standard model of cosmological infiation taking into account 

ff^ ' extrinsic curvature effects. We show that the modified model allows for a solution 

where the scaling function vanishes like a power law with exponent (\/3 + l)/4 ~ 
• i 0.683 at some initial time. 
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1 Introduction 



^ ■ A fundamental postulate in Einstein's theory of gravitation is that physics only depends 
- - on intrinsic geometric properties of curved space-time. Many approaches to gravity and 

cosmology assume that space-time is a brane, i.e., a submanifold embedded in a higher 
dimensional ambient space [T] (see |2] for an extensive list of commented references). 
In such a brane world scenario, Einstein's postulate cannot be assumed but rather is a 
property that should be derived from an underlying theory. One intriguing possibility is 
that this property holds true only in some approximation. In this case one might be able 
to predict physical effects that depend on extrinsic curvature. This would be interesting 
as a means to falsify theoretical proposals. It also might suggest experiments that give 
direct evidence for (or against) the existence of extra dimensions. In this paper we propose 
a possible cause for extrinsic curvature effects that, as we believe, is potentially relevant 
in brane world scenarios. We also discuss possible physical implications. Our proposal is 
motivated by well-established results on constrained quantum mechanics, as discussed in 
the next paragraph. 
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As proposed already by Schrodinger in 1926, it is natural to postulate that the quan- 
tum mechanical Hamiltonian for a free particle on a submanifold of three dimensional 
Euclidean space is proportional to the Laplace-Beltrami operator on this submanifold, 
and it thus only depends on intrinsic geometry; see Equation (31) in [3j. A more physical 
method is to derive this Hamiltonian as follows: use the well-established Hamiltonian 
in three dimensional space and restrict the particle to the manifold by a strong con- 
fining potential. As discovered by Jensen and Koppe [3] in this context the effective 
Hamiltonian thus obtained contains, in addition to the expected Laplace-Beltrami oper- 
ator term, an induced potential that depends on intrinsic- and extrinsic curvature; see 
[SI m El 13 [ini [m 1121 [IS] for generalizations and alternative derivations of this result. For 
example, the induced potential for a curve is attractive and proportional to the curvature. 
This explains why, as shown by Exner and Seba an electron on a wire has bound 
states when the wire has non-zero curvature]^ see also |15] for similar results in the con- 
text of electromagnetic waveguides. The existence of the induced potential was confirmed 
experimentally; see e.g. [TB] and references therein. 



Consider a real- valued Klein-Gordon field cj) on curved space-time embedded in higher 
dimensional Minkowski space. We propose to derive an effective equation of motion for 
(f) from a Klein-Gordon field on ambient space restricted to the space-time submanifold 
by a strong confining potential. As we show, this yields an effective description by the 
following modified Klein-Gordon equation on space-tim^ 



\g\-'/%\g\'/'9'"'d, + VU^)] + (moc)2)0(a;) = (1) 



with the induced potential Vind in Equation f ll4p below (we use Planck's constant h and 
the vacuum velocity of light c only here, to emphasize that the induced potential is a 
correction to the Laplace-Beltrami operator; in the rest of the paper we set ^ = c = 1). 
In fact, our result is more general: we use the Lagrangian formalism, and ([1]) corresponds 
to the Euler-Lagrange equations of the non-interacting part of the action only, but our 
result applies also to the case with interactions. It is worth noting that Vind is universal 
in the sense that it is independent of the details of the confining potential. We also 
investigate if and when Vmd could lead to measurable effects. For that we compute Vmd 
for various special cases. 

We emphasize that, for given space-time, the embedding in a higher dimensional 
Minkowski space is not unique, and the induced potential depends on the embedding. Our 
proposal can therefore only lead to a definite prediction in combination with a theory that 
fixes the embedding. We do not assume such a theory but instead use simple embeddings 
in flat ambient spaces with the lowest possible dimension. Fortunately, in many examples 
of physical interest (including the ones we discuss), there exists one such embedding which 
is preferred by its naturalness and simplicity; see [17] for a list of known examples. We 
thus believe that our approach is justified by "Ockham's razor". The investigation of 
more complicated embeddings (motivated by string theory etc.) is left to future work. 



similar result was found earlier by Marcus [5] in an effective quantum model for chemical reactions. 
^Note that the intrinsic geometry of a curve is trivial, and thus Schrodinger's postulated effective 
Hamiltonian cannot explain such bound states. 
■^Our notation is explained in the next section. 
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While the formula for the induced potential we obtain is identical with the obvious 
generalization of the well-known one in constrained quantum mechanics [H [3 |8] from 
Riemannian- to pseudo-Riemannian spaces, the physics is very different, and the deriva- 
tions of the induced potential in the literature do not apply to our case. We therefore 
present a different derivation for pseudo-Riemannian manifolds Ai"^'^ of arbitrary signa- 
tures {q,n) embedded in flat ambient spaces M^-'^+p (i.e., the number q > and n > 1 
of time- and space-like dimensions in M'^''^ is arbitrary, and so is the number p > 1 of 
the extra dimensions). To not further burden our notation, we restrict our derivation to 
the case g = 1 of main interest to us, but the generalization to arbitrary q is obvious. To 
be more specific: Derivations of the induced potential in constrained quantum mechanics 
are usually based on equations of motion; see [TH], Section LA for a more detailed discus- 
sion of the history and the different levels of mathematical rigor of these derivations. A 
key point in such quantum mechanical derivations is the probability interpretation of the 
quantum mechanical wave function. This provides an argument to use a scaling factor for 
the wave function in ambient space, and this scaling factor leads to the induced potential; 
see e.g. [TT], Equation (12) if for a lucid discussion of this point. However, Klein-Gordon 
fields do not have such a probability interpretation, and this argument therefore cannot 
be used. Instead, we start with the standard action for a Klein-Gordon field on fiat am- 
bient space and with a suitable strong confining potential; see Equation ([5]) below. By 
expanding this field in suitable modes and straightforward computations, we find that 
this action describes coupled Klein-Gordon fields on space-time; see Equation (ITT]) . Using 
a standard physics argument, we finally reduce the latter to an effective action for a single 
Klein-Gordon fields; see Equation (|T7j) . The above-mentioned scaling factor arises in this 
derivation for purely mathematical reasons. Our derivation suggests that such an effective 
action cannot be obtained in cases where the number of time-like directions in physical 
space-time and ambient space are different. 

In the rest of this paper we first introduce our notation and present our derivation 
of the induced potential for confined Klein-Gordon fields (Section [2]). Our examples are 
in Sections [3HS1 We first give the induced potential for the Schwarzschild metric and 
discuss its possible physical implications for black holes (Section [3]). We then consider 
the induced potential for the closed Robertson- Walker metric (Section Hj), and we use this 
to construct and study a brane-world model of cosmological inflation taking into account 
extrinsic curvature effects (Section E]). We end with conclusions in Sectional Some details 
of our computations are given in two appendices. 

We will mention basic results on general relativity and cosmology that are discussed 
in several textbook (including [151111]) without further reference. 

2 General result 

We assume that space-time Ai = A^^'*^ is a (n -|- l)-dimensional Lorentzian manifold 
embedded in {n + p + l)-dimensional Minkowski space M^'^+p. 

Notation: We use capital latin letters M, N, . . . for indices running over 0, 1, 2, ... , n+p, 
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and denote as Z = {Z^ , , Z"^ , . . . , Z"''^^) inertial coordinates in R^'^+p, i.e., the line 
element in these coordinates is (the following defines our sign conventions) 

ds^ = 7]MNdZ^'dZ'' = {dZy - {dZ^f (2) 

Greek letters /i, z/, A, a, . . . and lower case latin letters i,j,k... are used for indices running 
over 0, 1, 2, . . . , n and n + l,n + 2, . . . ,n + p, respectively. We write bold face letters for 
vectors in R-'^-'^+p and a dot for the scalar product of such vectors, i.e., v • w = rjMN'v'^^'w^ . 
We assume that space-time Ai can be (locally) parametrized by a function Z = f (x) with 
X = (x°, . . . , x") in an open subset of R^'". The metric tensor of M. is then g^^, = ■ t^, 
with tangent vectors t^{x) = 5^f(x)0 Using y = {y'^~^^, . . . ,y"'~^^) € R^ we construct a 
local coordinate system of an open neighborhood oi A4 G R^'""'"^ by setting 

~f{x,y) = f{x)+fn,{x) (3) 

with p linearly independent space-like vectors rij orthogonal to all tangent vectors and 
with constant scalar products, i.e., 

t^(x) ■ ni{x) = 0, ni{x) ■ nj{x) = -h^ (4) 

for all X, fi and i, with constant hij defined by this equation. We can further restrict 
the vectors rij by the additional condition hij = 6ij, without loss of generality. However, 
we sometimes find it convenient to leave hij general in our equations. We assume that 
there exists e > such that the coordinate system (x, y) is one-to-one for all y such that 
\y\ = \/hij]/y^ < e (as discussed in the beginning of Appendix lA.H we believe that this 
is a minor restriction). We denote the metric tensor on R^'^+p in the coordinates {x,y) 
by GMN{x,y). Here and in the following we write \h\, \g\ and \G\ for the absolute values 
of the determinants of (hij), (g^u) and (Gmn), respectively. We define G^^^ as usual: 
G^^Gln = 5%, and similarly for gf"" and h'^ . 

We consider a Klein-Gordon field $ = $(Z) on ambient space R^'"+^' with the usual 
dynamics but confined to the region close to the submanifold by a potential Vcon/(Z), 
i.e., the (free part of) the action is 

^0 = ^ / rf"+^+^z(r/^^(aM$) - (mLe + V^onf)^') (5) 

with the "bare mass parameter" m^^^g e R. We assume that the potential Vconf is strongly 
confining to M, i.e., in the coordinates {x, y) it has the form 

V,onAx,y) = Kon/(f(x,2/)) = ^Viy/e) + 0{\y\) (6) 

with e > a convenient scaling parameter assumed to be small |f| we allow for a possible 
correction term 0(|y|) depending on x but vanishing at least linearly with \y\ as \y\ — > 0. 

''We write Om, 9^ and di short for -q§w, gfrr and respectively. We always assume implicitly that 
functions we introduce are differentiable up to the degrees needed. 

^Introducing this scaling parameter is a useful mathematical trick allowing to cleanly separate different 
energy scales; see e.g. [13] for a lucid discussion of this point. 
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We assume that V{y) is such that the eigenvalue equatioi 

(y) (7) 

has a unique solution xo corresponding to the smallest possible eigenvalue /xq and such 
that / d^y \h\^^'^\y\'^\xo{y)\'^ is finite for all = 0, 1, 2, . . ., and fia — fJ'O > ^ for all a ^ 0; 
here and in the following, indices a, P . . . denote quantum numbers labeling the solutions 
of the eigenvalue equation in ([7]). Below we refer to such a potential V{y) as suitable. We 
choose the eigenf unctions Xa to be real-valued and normalized such that 



J d^y \h\'/\Mxp{y) 



To be specific we mention one example for a suitable potential allowing for simple com- 
putations of all solutions of ([7]): hij = 6ij and 



n+p 



Viy)= hm vo(l- n eii-\y^)) (9) 

j=n+l 

with i > and the Heaviside function 6 ("infinite box potential") but, as discussed in 
Appendix \A.1\ the same result is obtained for a large class of potentials. In particular, V 
can be bounded and such that the eigenvalue equation in ([7]) also has scattering solutions 
(in this case the symbols and Sap have to be partly interpreted as integral and Dirac 
delta). 

One can expand $(x, y) = ^{i{x,y)) in the eigenf unctions e^^^'^Xaiy /^) and thus 
rewrite the action in (|5]) as an action of an infinite number of fields 0a (x) on space-time 
A4. As explained in Appendix IA.2t a key point in this computation is that, to correct a 
mismatch of Jacobian determinants, one has to include a scaling factor (|/i| l^fl/IGI)^/^ in 
this expansion as follows, 

Hx, y) = ( E <t^ai^)^'"'x.{y/e). (10) 

By a straightforward computation we obtain the following. 

Result: For suitable confining potentials V and sufficiently small e > 0, the action in ^ 
equals 

with 

c,.p = V / ^'y \h\'^'y'xp{y)daa{y) (13) 



^To be more precise, we assume that V{y) is such that the fohowing is an eigenvalue equation of a 
self-adjoint Schrodinger operator on the Hilbert space of square integrable function on MP. 
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and 

V^nd = \h'^ - 2a.A'^«,/) (14) 

where Oiii^{x) and A-J'lx) are defined by the following equation, 

d^Ui = -a,/tA - A./nfc. (15) 

(See Appendix lAl for a derivation of this result.) 

We note in passing that the action in f lTT]) - f|T4|) has a gauge theory structure discussed 
in [HI [11], for example, but this is not used in the present paper. 

Thus, up to terms 0(e) vanishing in the limit e — )■ 0+, the confined Klein-Gordon 
action on ambient space in ([5]) is equivalent to an action of Klein-Gordon fields (pa on Ai 
with effective masses 

ml = mlare + f^a/e^ (16) 

and additional potential- and derivative terms. It is natural to fix the renormalized mass 
ml and choose m^^^g = — /ig/e^. The other mass parameters m^_^g = + (/Xq, — /io)/e^ 
are then much larger than ml. Standard physics arguments suggest that, for sufficiently 
small values of e, only the Klein-Gordon field 0o with the "small" mass parameter m^ 
is relevant for the low energy physics properties of the model. It therefore is a good 
approximation to simplify the model by replacing the action in ffTTj) with 

^eff = \ I d'^'-'x \g\'/' {g''''{d,<P) (9.0) - (ml + F,„,)02) (17) 

and (p = 00- This is our result described in the introduction. It is important that it is 
independent of the details of the confining mechanism: changing V can be compensated 
by a change of the bare mass parameter m^^^g and thus is irrelevant. 

As already mentioned, the generalization of this result to interacting Klein-Gordon 
fields is straightforward. In particular, changing the action in ([5]) by a ^^-interaction 
term, 

S = So-^ [ rf"+^+lZ<l>^ (18) 
leads to the following change of the effective action in f[T71) . 

^e// = ^S/-T / (19) 
^ Jm 

with the renormalized interaction strength 

Ao = A,,,ee-P / d^y\h\'/\oiy)'. (20) 

Again it is natural to fix Aq > and adapt Xbare accordingly, i.e., the effective action is 
independent of the details of the confining mechanism also in the presence of interactions. 
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3 The Schwarzschild black hole 



The Schwarzschild metric is a spherically symmetric solution of Einstein's equations de- 
scribing a black hole with mass M in an otherwise empty universe. It is given by the line 
element 

ds^ = (l - ^) de - (l - ^) ^ dr^ - r\de^ + sin\e)dip^) (21) 

where = 2GM is the Schwarzschild radius and x = {t, r, 6, ip) the Schwarzschild coor- 
dinates, as usual. 

As proved by Kasner [20] , the Minkowski space-time of lowest dimension allowing an 
embedding of Schwarzschild space-time is M^'^. One such embedding found by Fronsdal 
is given by 

2rs(l - r^/r)^/^ sinh(t/[2rs]) (r > r^) 
2r,(r,/r - 1)^/2 cosh(t/[2r,]) (0 < r < rj 

2'rs(l - rs/ry/"^ cosh(t/[2r^,]) (r > r^) 
2r,(r,/r - 1)^/2 sinh(t/[2r,]) (0 < 
= g{r) 

= r sin(^) cos{ip) 
Z"^ = r sin(^) sin(y9) 



Z^ 



(22) 



Z^ = r cos{e) 



where 



^ ^ ~ dr \l r3 ■ ^ ^ 

This embedding is natural in that no other complete embedding of the Schwarzschild 
solution in the six-dimensional Minkowski space-time exists [22]. 

By straightforward computations we find the induced potential (see Appendix IB . 1 1 for 
details) 

(f3 + f2+f + 9)(f2+ l)(f + l) . , 

Vindir) = ^^A,J ^ T ^ r = r/rs (24) 

16r^r*(r"' + r + 1) 

which is non-singular in the whole region < r < oo. Considering that the details of the 
embedding are different for r > and r < r^, it is remarkable that the induced potential 
is continuous at r = r^. Note that Vindi'f') is strictly monotone and has the following 
asymptotic behavior, 

' -9/(16r2f^) (l + f/9 + 0(f2) j (r ^ 0+) 

Vind{r) = l-l/r',{l-3{f-l) + 0{{f-lf)^ (r ^ r,) . (25) 
-l/(4r,)2(^l + l/f + 0(l/f2)) (r^oo) 

Moreover, it approaches its limiting value — l/(4rs)2 rather rapidly and, for distances r 
larger than lOOr^, Vind{f) differs from its limiting value by less than 1%. Note that the 



7 



magnitude of this limiting value corresponds to a boson mass, in SI units, 

_n_ ^ 3.31xlO^°kg GeV 

4cr, ~ M c2 ■ ^ ' 

We conclude this section with a short discussion of the possible physical relevance of 
our result. It is worth noting that the boson-mass equivalent of the magnitude of the 
induced potential in fl26l) is proportional to the inverse of the black hole mass M. More- 
over, it seems that only very light black holes, with masses 10^" kg or smaller, could give 
rise to induced potentials that have measurable effects, and the spatial variation of this 
potential is significant only very close to such a black hole. Thus the natural candidates 
to search for such effects are primordial black holes; see e.g. [23] and references therein. It 
is remarkable that, very far from a black hole, the induced potential renormalizes boson 
masses by a negative constant. If many primordial black holes exist in the universe their 
cumulative mass renormalization effect could be quite large. 



4 Robertson- Walker space-time 

The Robertson- Walker space-time describes the evolution of a homogenous isotropic uni- 
verse. Its metric depends on a real parameter K whose sign determines if the universe 
is open {K > 0), flat {K = 0), or closed {K < 0). All cases can be described by a line 
element of the form 

= - a{tf {dr^ + S{rf{de'' + sin^{e)d^^)) (27) 

in coordinates x = {t, r, 9, ip), with t the cosmological time, a{t) the so-called scale factor, 
and 

{sm{y/Kr)/VK {K > 0) 

r {K = 0). (28) 

sinh(v^r)/v^ {K<0) 

This space-time can be naturalljfl embedded in M}'"^ as follows (see [17], B.5.), 





r b{t)/VK 


{K > 0) 




(l/2)(rVro + ro)a(t) + 5(t) 


{K = 0) 




[a(t)C(r) 


{K < 0) 




ra(t)C(r) 


{K > 0) 




(l/2)(rVro-ro)a(t) + 5(t) 


{K = 0) 




[Kt)/V\K\ 


{K < 0) 


Z' = a 


,(t)S{r) sm{9) cos{{p) 




Z^ = a 


,(t)S'(r) sin(^^) sm{(f) 




Z' = a 


it)S{r) cos{e) 





(29) 



^We believe that there is no other embedding in five dimensional Minkowski space-time, but we are 
not aware of a proof of this in the literature. 
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with functions b and B defined by the following equations, 



b{t) = ^K + a{ty {K 0) (30) 

and 



Bit) = TT (31) 



{b{t) = db{t)/dt etc.), and 



C(r) ^ > 0) . ^32) 

^ ' ^ cosh(v^r)/v^ (ir < 0) ' ^ ' 



ro > is an arbitrary parameter. 



By straightforward computations we find the induced potential (see Appendix IB. 21 for 
details) 

- 6- + 3 5 — — .„ 

4 V a a? K + a? 



V.n,= -{6- + ^—;^-j—-,] (33) 



for arbitrary K. Considering that the details of the embedding are different in the three 
cases K > 0, K = 0, and < 0, it is remarkable that the induced potential has such a 
simple form, and, in particular, that it is continuous at K = 0. 

We conclude this section with a preliminary discussion of the possible relevance of 
this induced potential in cosmology. To simplify some formulas we only consider the case 
K = 0. Inserting H = a/a we find by simple computations 



V^ind = -^(g' + 6g-3) (34) 

with q = aa/d^ = —{H + H^) the usual deceleration parameter. Thus Vind > if the 
deceleration parameter is in the range —6.46(4) < q < 0.464(1) and < otherwise, and 
Vind has typically the same order of magnitude as H^. Since the value of the Hubble 
constant Hq today corresponds to a boson mass of approximately 10~^^ eV, this suggests 
to us that the induced potential is negligibly small in the universe at present and far back 
in time. To see if the induced potential could have an effect at early times we compute 
Vind for a universe with a scale factor growing as a(t) = {t/t^Y some exponent x > 
and some constant t*, with t the time after the big bang at t = 0. We find 

- 4x - 1 

Vind = , (35) 

i.e., as t — )■ 0+, the induced potential diverges towards — oo for < x < (1 + -\/3)/4 and 
towards +oo for x > (1 + a/3)/4. This suggest that the induced potential could have 
had an important effect in the early universe. In the next section we propose and study 
a self-contained model for the evolution of the early universe, taking into account the 
induced potential. As we will see, this model predicts that, as t — )■ O"*", the scaling factor 
can vanish like a{t) — )■ {t/t^Y with x = (1 + "\/3)/4 the critical value of the exponent, and 
such behavior is impossible without the induced potential. 
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5 Extrinsic curvature effects in the early universe 



Models of the early universe often assume a Robertson- Walker metric and matter repre- 
sented by an isotropic Klein-Gordon field; see e.g. [211 [25] for reviews or [19j for a recent 
textbook on this topic. In this section we propose a generalization of such a model taking 
into account extrinsic curvature effects. We also present results on the solution of this 
model and, to put them in perspective, compare them with results for the corresponding 
standard model where the external curvature effects are ignored. 

We obtain this model by adding our induced potential term to the action consisting 
of the usual Einstein-Hilbert term, cosmological term, and Klein-Gordon term with 0^- 
interaction |19]1^ 

S=lld'x \g\'/'(^MUR - 2A) + g^''{d,<P){d.<P) - {m' + xV,nd)<p' - ^<P') (36) 

with the Ricci scalar R and the induced potential Vind in f|T4l) . Our model parameters 
are the reduced Planck mass Mpi = I/a/SttG, the cosmological constant A, and the 
renormalized mass m and coupling constant A of the Klein-Gordon field, respectively. 
To avoid writing similar formulas twice we use a parameter x which is either (for the 
standard model case) or 1 (for the extended model case with extrinsic curvature effects 
included) . 



5.1 Model details 

We are interested in a homogeneous Klein-Gordon field (f) = (f){t) (depending only on 
cosmological time t) on Robertson- Walker space-time with scaling factor a{t). The stan- 
dard method to derive differential equations determining the time evolution of a and 
in the case x = is as follows: restrict the Euler-Lagrange equations obtained from the 
action in (!36|) to the Robertson- Walker metric in (127]) and Klein-Gordon fields = 0(t). 
Unfortunately we cannot use this method in the case x = 1- this would require a formula 
describing how the induced potential changes with arbitrary variations of the metric, but 
this we do not have (since the induced potential depends on the embedding, and we do 
not know how to naturally change the embedding with the metric in general). 

We thus use an alternative method which, in the standard case x = 0, leads to the 
same result as the method just described: We insert the metric in (127]) and = 0(t) 
into the action in f l36p . This yields S = const j dt L (we can ignore the multiplicative 
constant) with the Lagrangian 



L = M''pi{-3aa' - Aa') + — i - (m^ + ^ Q- + 3^ - ^ )0^ - -0^ 1 (37) 



2 \ " V 4 



• 2 ■■ 2 

a a a 



a a? a? 



X 



where we used fl55]) . setting K = for simplicity (we dropped total derivative terms; 
a = da/dt etc). Note that, in the standard case x = 0, the Lagrangian in (137]) is of the 



^In this section we write A and m short for Aq and toq, respectively. Moreover, the meanings of the 
symbols x, y, and 0„ are different from the ones in Section [2] 
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following usual type in dynamical systems: L = L{a,a,(j),(f)), corresponding to second 
order differential equations familiar from mechanics. However, in the case x = 1? the 
Lagrangian also depends on a, and this leads to a 4-th order time evolution equation for 
a; see e.g. for some general background on such higher-derivative Lagrangian systems. 
We note that, even though the results described below are well-known in the case x = 0? 
our method seems different from the ones used in the literature. 



By straightforward computations we obtain from fl57|) the following Euler-Lagrange 
equations, 

if2 



and 



+ 3H(j) + + x 



H + 2H' 



+ = 



2ii + 3fi2 _ A + J_ f- - —0' - --^^ 



M2 ^ V 2 2 



A 

T 



+ 



+ 



X 



H 



HH 5H SH"^ 



2H AH^ 2 
if2 



+ 



H 



H 2H H 



3H^ 12 8H^ 12H^ 



H 



(3^ 



(39) 



+ m + 



with H = a/a. It is interesting to note that only H and its derivatives appear in these 
equations, and thus the order of the Euler-Lagrange equation corresponding to a is effec- 
tively reduced by one. It is important to note that the Lagrangian in ( 1371) is invariant 
under time translations, and therefore the time evolution equations in fl38|) and fl39|) allow 
for a conservation law given by [26] 



da dt \ da 



.. dL -dL ^ 

flTTTT + - L 

oa d(h 



(40) 



(using the Lagrange equations one easily checks that dl/dt = 0). By straightforward 
computations we find that, up to a factor a{t)^, this conservation law can be expressed in 
terms of H and its derivatives as follows: 



3a3M|,; 



ml, 



m 



2 J,2 



A 
2' 



+ 



X 



6H 3 . 



+ 



IH^ 1 H 1 „ 1 . 

H + -H 

8H^ 12 H 6 4 . 



(41) 



In the standard case x = 0? the two Friedmann equations are equivalent to 
and X = 0: the dynamical system defined by the Lagrangian in flTTI) allows for many 
more solutions than Einstein's equations obtained by varying the action in f l5B]) . but the 
additional solutions are eliminated by imposing that the constraint that the value of the 
conservation law in (HOj) is zero. We assume that this is true also in the extended case 
X = 1- We emphasis that this is a hypothesis but, as we believe, a plausible one. To 
check it one should generalize the method described in the previous paragraph to the case 
X = 1, but this we leave for future work. 
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We thus propose the equations in fl38|) . fl39|) . and X = with X in fHTj) and for x = 1, as 
extension of a standard cosmological model by extrinsic curvature effects. We are mainly 
interested in this model short after the big bang. We thus make the following ansatz for 
solutions, 

0(t) = (00 + 0lt + 02t' + 03^' + ■ ■ ■ ) 

X (42) 
E{t) = - + Ho + H^t + H^f + Hst^ + --- . 

We call solutions with x = generic, and the others scaling solutions. The reason for these 
names is as follows: since the equation of motions are invariant under time translations, 
one can replace in (H2|) t hj t — to, to arbitrary, and get another solution. Generically, the 
fields a{t) and 0(t) at their derivatives at some time to are finite, and the generic solution 
gives a series representing these fields in the vicinity of t = to in terms of appropriate 
initial conditions at t = to- For scaling solutions, on the other hand, the time t = (say) 
is special: it is the time when the scaling factor vanishes like (t/t^)^ for some > 0, and, 
if y > 0, the boson field diverges. 



5.2 Solutions 

We now discuss solutions of , (15^ , and X = with X in (HTi) , for the two cases x = 
and X = 1 (some details on how we obtained these solutions are given in Appendix IB. 30 . 

For X = 0, we found the following generic solution, 

(Pit) = 00 + (pit - ^{SHo(Pi + K + A0^]0o)t' + 0(t=^) 

Hit) =Ho~ ^t + ^iSHoh + K + A0^]0o)t^ + Oit') (43) 



= y^A/3 + (0?/2 + mm/2 + A04/4)/(3M' 



2 1 
PI) 



(the coefficients of the 0(t^)-terms are given in Appendix IB. 31 (159!) ). This solution de- 
pends on two free parameters 0o = 0(0) and 0i = 0(0), and this is the maximum possible 
number of free parameters: for x = 0, (!38l) and (!39|) are second- and first order differ- 
ential equations, allowing for two- and one integration constants, respectively, and the 
constraint X = reduces the number of free parameters by one. 

In the standard case x = 0, we did not find any solutions as in with x 7^ 0. 
However, in the extended case x = 1; we found two such solutions: one with x = 3/4, 
y = I given by 



(44) 
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and another one with X = (1 + V3)/4 = 0.683(0), y = given bf] 

0(t) =Mp/l.50(2)- [o.228(0) + 0.402(0) AM|; + 0.0883(5) aI t2 + 0(t 



i7(t) 



0.683(0) 



- 0.0142(5) + 0.120{0)XMpi - 0.0690(4) A]t + 0{r). 



(45) 



Note that the solution in (l44l) exists only if A > 0, whereas the solution in (l45l) exists even 
if A = A = 0. Thus the latter solution is more robust and, as we believe, more interesting. 
We therefore give formulas for the coefficients of the 0(t^)-terms in (H5ll in Appendix IB. 3 

For X = 1 we found the following generic solution, 



(f){t) =00 + 0it - 2 ( 3^001 + 00 



tt2 

+ X<Pl + Hi + 2Hl - ^ 



+ O(i') 



H{t) =Ho + H^t+{ + ^02 + _ 



^0 + 



3Hf 
AH, 



(46) 



+ - H,)^ + (0? + 2AM2^ - 6//2m|,)^ + 0{t 



which depends on four free parameters: 0o ^ 0, 0i as before and, in addition, Hq = H{0) 
and Hi = H{0). Note that this is the maximum number of free parameters (since ( l39l) 
now is third order, we have two more free parameters). 

We note the static solutions, i.e., 0(t) = 0o and H(t) = Hq independent of t, which 
are interesting special cases of the generic solutions above and exist if < 0, A > 
("Mexican hat potential") and A > m''/(4AM|,J. They are given bj0 



Ha 



4AA - mVM 



PI 



12A 



(x = o) 



(47) 



and 



/2(-3m2 - 2A) 
6 A + mVM|,; ' 



Hn 



4AA -mVM 



PI 



2(6A + mVM' 



piJ 



(X = l) 



(45 



in the standard- and extended cases, respectively. Note that, for x = 1, we have the 
additional conditions A < — 3m^/3 and A > — m^/(6Mp;) for these solutions to exist. 
Moreover, for A > and in both cases x = s-iid X = 1) there exists also the static 
solution with 0o = and Hq = 



5.3 Discussion 



As already mentioned, the generic solutions allows us to compute a{t) and 0(t) in some 
interval around an arbitrary initial time t = to from its initial conditions at t = to- In 

^We write 0.683(0) for a numerical value 0.6830 ± 0.0001. 

^°Note that for every solution (j){t), Hit) we give, —(pit), H{t) is also a solution. 



13 



the standard case, they allow us to understand the importance of the so-called " slow roll" 
condition for inflation in this model pi)J in the following way: the equations of motion 
of the boson field (f){t) are identical with Newton's equations for a particle moving in one 
dimension under the influence of a friction term and a potential V(0) = m^0^/2 + A(/)^/4: 
4){t) = -3H{t)^{t) - V'((/)(t)), and implies H{t) = -0(t)V(2M|,;) and H{t) = 
—6H{t)H{t) + dV {4>{t)) / dt (these equations can also be obtained directly from the Fried- 
mann equations). This suggests that, starting from arbitrary initial conditions 0(0) and 
0(0) at some time t = 0, (f){t) will increase with decreasing 0(t) until (f){t{)= at some 
time t = ti, and, at this time, H{ti) = H{ti) = 0. In the vicinity of t = 



0(t) = 00 - ^V'(0o)(t - + 0((t - ti)2) 



(49) 



with 00 = 0(^1 ) and Hq = a/A/3 + V(0o)/(3Mp;), i.e., the time evolution of the scaling 
factor can be well approximated by the exponential law a{t) ?a a{ti) exp[H{ti){t — h)), 
up to corrections that are negligible as long as 0(t)^ <^ 24:MpiH(ti)/\t — ti\. We computed 
and plotted the approximation of the solution obtained by extending this series in ( H9l) to 
higher orders and ignoring the 0{{t — ti)'^)-terms, for different values of n, up to n = 8. 
The result suggests that these series have a finite radius of convergence which, however, 
is often smaller than the full time interval where inflation occurs. 

We now discuss the extended case x = 1 how it differs from the standard case 
X = 0- Obviously, there are two essential differences: firstly, the generic solution in fH6!) 
depends on four (rather than two for x = 0) initial conditions, and secondly, there exist 
scaling solutions (which do not exist for x = 0). In the next two paragraph we discuss 
these two differences and possible implications. 

The generic solution in ( H6|) depending on four initial conditions is a consequence 
of the the Lagrangian depending also on d{t) As illustrated in a simple toy model in 
Appendix IB. 41 even if such higher-order term in the Lagrangian is very small, it can lead 
to a much richer qualitative behavior of the system. Thus, at first sight, it seems that the 
model with x = 1 is less predictive than the standard model with x = 0- However, our toy 
model also suggests that it is possible to restrict the initial conditions of the model with 
X = 1 so as to obtain solutions that are similar to the ones for x = 0. One natural way 
to impose such a restriction is suggested by the generic solution in fll^ : obviously, this 
solution is not well-defined if 0o = 0. Using translation invariance in time, we conclude 
that a generic solution such that 0(to) = at some time t = to can exist only if we impose 
some restrictions of the free parameters. One can check that these restrictions on initial 
conditions at time t = to are 

^2 _ A , 0? rr _ _ J0|(0i+2AM^ 



3 6M|,;' 3M|,j(30? + 4AM2 



PU 
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We checked that the foUowing result holds true for arbitrary differentiable boson potentials V ((/)). 
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with Ho = H{to) etc. (see Appendix IB.3I for derivation). This leads to 

(Pit) = - to) - l<PiHoit - to)' + out - to)''), H{t) = Ho + H,{t- to) 

, 4Hi(Pl{19(Pi + 360f M^,A + 32(AM|,)^) ,.2^n((, . ^^^^ 
+ {3(Pl + AAMj,,)^{5(Pl + 4AMl,,)) ' ''^ +O{{t-to)) 



with Ho and Hi as in ( 1501) . We thus propose the following hypothesis: Only solutions that 
remain well-defined when (p{t) as t to, for some time to, can describe the evolution 
of the universe. This reduces the number of free parameters to two, which is equal to the 
number of free parameters in the standard case x = 0: the time t = to where the boson 
field (p{t) vanishes, and the time derivative (pi = (p{to) of the boson field at that time. We 
believe that it is these restricted solutions of the extended case x = 1 that behave, at 
larger times, similar to the ones on the standard case x = 0. We note that there might 
be also solutions where (p{t) > for all times, but these we do not know how to restrict 
in a natural way. 

We believe that the existence of scaling solutions in the extended case x = which do 
not have any analogue for x = O5 could make the former model more predictive than the 
latter in the following way: for non-linear dynamical systems, scaling solutions often are 
attractors that capture the qualitative behavior of a large class of solutions (see e.g. [27]). 
As mentioned already above, the scaling solution in (H5l) is more robust and thus, as we 
believe, more interesting than the one in (jS]). We thus believe that it would be interesting 
to explore the validity of the following hj^othesis: Solutions that can describe the universe 
immediately after the big bang approach the scaling solution in ( f73P as t — t- 0+ . Combining 
this hypothesis with the one above we get the following picture: at some time to > 0, 
(p{to) = 0, and in some interval the time evolution is described by f l50l) - fl^Tl) . The free 
parameters (pi and to of this solution are then to be fixed such that, as t — )■ 0"*", this 
solution approaches the one in fH5|) . For t > to, the boson field (p(t) grows with decreasing 
(p{t) until (p{t) = at some time t = ti, and close to this time the solution is, as we 
expect, similar to the one in (H9|) . If or not this picture holds true can be determined by 
a numerical solution, but this we leave to future work. We expect that the solution close 
to the time t = ti describes inflation, similarly as discussed after (H^ . One important 
question is if such a solution can describe inflation (solving the horizon problem etc. jT9] ) 
even in some time interval short after t = or not, i.e., if such a solution has one or 
two periods of inflation. At first sight it seems the answer to this question is negative: 
a{t) — )■ {t/t^)^ as t — )■ 0+ with x = 0.683(0), and previous work on power-law inflation 
[2Sl[2n] suggests that inflation requires a; > 1. However, in our scaling solution, a{t) is well 
approximated by {t/t^)^ only in a small time interval, and for larger times the behavior 
is more complicated. We thus believe that this question is open. 



6 Conclusions 

In this paper we proposed that the induced potential, which is well-established in con- 
strained quantum mechanics, could also be relevant in brane-world scenarios where space- 
time is embedded in a higher dimensional ambient space. We showed that, by assuming 
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that the propagating degrees of freedom are restricted to space-time by a strong confining 
potential, one finds that the Klein-Gordon equation on space-time is indeed modified by 
a induced potential term which does not depend of the details of the confining poten- 
tial. While the formula for the induced potential we obtain is the same as in constrained 
quantum mechanics, the physical arguments to derive this potential are different. 

One striking feature of the induced potential is that it does not only depend on 
intrinsic- but also on extrinsic geometric properties of the embedded space-time. Thus, 
if the induced potential has observable consequences, it offers the intriguing possibility 
to test brane-world scenarios experimentally. As examples we computed the induced po- 
tential for Schwarzschild- and Robert son- Walker space-times. Our results suggest that, 
while the induced potential is usually negligibly small, it might be relevant in extreme 
situations like the early universe or in regions close to primordial black holes. We also 
proposed and studied a model for cosmological inflation with the effect of the induced 
potential included. At first sight this generalized model seems to be less predictive (since 
the generic solution of the equations of motion depend of four initial conditions, rather 
than two in the standard case). However, we found a natural condition reducing the 
number of initial conditions to two. Moreover, the extended model is different from the 
standard one also in that it allows for scaling solutions that do not depend on any free 
parameter. These scaling solution describe an expanding universe with a size approaching 
zero at some initial time. Our results suggest that, in the extended model, one only has 
to postulate a big bang, and the time evolutions of the universe is fixed, without much 
freedom to vary initial conditions. Anyway, it seems worthwhile to study this model 
further. 

For simplicity we restricted ourselves in this paper to the case where the ambient 
space is fiat. It would be interesting to generalize our results in Section [2] to cases where 
the ambient space is curved. We expect that the induced potential thus obtained is 
the obvious generalization of one known in constrained quantum mechanics; see e.g. |9j. 
Equation (3.37). 

A Derivation of the induced potential 

In this appendix we give details on how to derive the result stated in Section |2i 
A.l Technicalities 

For compact manifolds Ai embedded in some ambient space, the existence of the coordi- 
nate system (x, y) as described in the beginning of Section [2]is guaranteed by the so-called 
Tubular Neigborhood Theorem. For this reason we believe that the assumption that such 
a coordinate system exists is not very restrictive. 

In the main text we argue that the details of the confining potential Kon/ are not 
important, and all that is needed is that this potential is such that the Schrodinger 
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equation in ([7]) has a non-degenerate groundstate solution, obeying rather mild decay 
conditions, and with a finite energy gap to the first excited state. While this is plausible 
from a physics point of view, we emphasize that our derivation of the induced potential 
below can be easily promoted to a mathematical proof only in the special case where the 
confining potential is (essentially) the box potential in (the careful reader will find 
the spots where we implicitly make this assumption). It would be interesting to find a 
mathematical proof that would apply to a much larger class of potentials. One technical 
challenge is to find arguments that avoid the use of local coordinates. 

A. 2 Derivation 

Taking the dot product of (fT5|) with t^, and recalling (jl]) we obtain 

(yk^lu = o^kfM^gxu = rifc ■ df,dj = Uku^ (52) 

which is useful in computations of the induced potential in examples. Similarly, by taking 
the dot product of (fTSll with iij, 

-^i^ij = ^kj = % ■ d^xii = —Aj^i. (53) 



We compute the action in (fTTj) by changing to the coordinates Z = (x, y) in M^'^+p 
defined in (|3]), i.e., Z'^ = x'^ for /i = 0, . . . , n and = for j = n + 1, . . . ,n + p. We 
obtain 

^0 = ^ / "^^^'"^ j ^"?/|^^r^'(^^'^'^(^Af$)(5^$)-(mLe + V;on/)$') (54) 

with ^{x,y) = <l>(f(x,?/)), Bm = q§m, and Gmn = (<9Mf) ■ (<97vf). Moreover, 

\ _ ( Gfj^u Gfj^j \ f 7^iy — G ^kh^^Gip Gfj^j \ 

with G^j = Gjf, = y^Ak^j and 

ipiu = 9fw - '^^y^otk^.y + y^y^otf^^^ai^)^. 

This allows us to compute the following, 

|7| = |det(7,..)| = -2y'=a,/ + yV(2a.AV'^-«A.A'^«//) +0(|i/|^)). (55) 



Due to its special form it is possible to compute the determinant and the inverse of 
the matrix {Gmn) in a simple way as follows, |G| = \'~i\\h\ and 

~ ' G^'^Gxkh^^ -h'^ + GkxGi„G^"h^'W^ 
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with (G^"^) the matrix inverse to (7^;/). Thus 

G^'^ = g^'' + 0{\y\) 

G^^ = = y V'^fcA' + Oi\y\') (56) 
G^J = _/,^. + + 0{\y\'). 

It is straightforward to compute the higher order terms in y^ but, as will be seen below, 
they are not needed. 

As mentioned in the introduction, it is important to rescale the Klein-Gordon fields 
so as to change the Jacobian IGI^^^ in (|5^ to We thus introduce = 171^/^$ 



So = lJ d"^'x\g\'/'J dPy\h\'/'(^G''''{dM<P){dN<p) - (mL, + Ko„/ + M*^') (57) 



with 


\l\ 


into 


m 




~ 2 


with 






4 



G 



MN 



^(^M|7l)(^iv|7l) - {dM~dN\l\)\ - \9\-'^\dM\l\)~dN{\9\"'G'''') ). 



Using fl55l) and fl56|) we find 

Vind{x,y) = Vind{x) + 0{\y\) (58) 

with Vindix) in dUD. 



We now can insert the expansion 0(x, y) = 4>a{x)'Xa{y) with Xa{y) = |e|~^^^Xa(l//e) 
and Xaiy) the eigenfunctions defined by (171) and (IHl). Note that the scaling is such that 
the Xa{y) are orthonormal. In the following we compute the different contributions to the 
integrand in (1571) . In our computations we encounter integrals 

d^y \h\^/Y' ■ ■■y'-{d^xM){d]xp{y)) = 0{e--^-') 

for a, 6 = 0, 1. Thus, for fixed a + 6, we only need to take into account the terms with 
n < a + h. This explains our truncations in fl56ll . 



Using (156]) we obtain 

d^y\h\^/^G^''{d^(l>){d,(l)) = Y,9^''idMidu(l^c.) + 0{t) 



and 

2 j dPy\h\^'^G^'i{d^ci)){d,<P) = 2Y,9^\dMCxo.p(t^p + 0{t) 
with Cxajs in ( !T3l) . The following integral is computed using ([7]), 

I d^y W^{-h^^{d4m,) - {mire + Kon/)0') = " J^^.e + ^^Je')€ + 0{e) 
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with the 0(e)-term coming from the 0(|?/|)-correction in (j6]). Using ( 156|) and inserting 

i = E,Id^y'H^'x,iy)x,iy') yields 

f (Fy \h\^'\G'^ + h'^){d.MM) = Yl 9''"C,,^C,,^(Po.^p + 0(e) 

with C^ai3 in ( fT3l) . Finally, recalling ( l58l) . 

- /" rf'^y|/ir/2V^,„,0' = -V^.„rf 5^0^ + 0(e). 

Collecting all integrals above we obtain the result in f llip - (IT^ . 

B Computation details 
B.l Schwarzschild black hole 

It is convenient to introduce normalized vectors that allows to write the embedding in 
(ESD as 

Z = f(x) = 2rs\l - rs/r\^^'^e^ + g{r)e2 + re^ 

for X = (t, r, 9, ip), i.e., 

^ ^ f (sinh(t/[2r,]), cosh(t/[2r,]), 0, 0, 0, O) (r > r,) 

[ (cosh(t/[2r,]), sinh(t/[2r,]), 0, 0, 0, O) (0 < r < r,) 

62 = (0,0,1,0,0,0) 

= (0, 0, 0, sin(6') cos(v9), sin(6') sin(y9), cos(6')). 

We find the following vectors normal to all tangent vectors and such that hjf^ = 6jk for 
J, k = 4,5, 

n4 = -g'{r)\r/r, - l\^^'^e^ + {rjrf''^e2 

n5 = {rs/rf''^\l - rs/r\^''^e^, + g'{r){l - rjr){r/rsf''^e2 - [rjrf'^er 

with g'{r) in ( l23l) . By straightforward computations we find the following non-zero com- 
ponents of a^^ , 

"40° = (^^ + + r2)V2/(2rr,), a^^^ = 3rj{r^ + rr^ + rl)"^'^ / {2r'^) 
(unless indicated otherwise, formulas hold true for < r < oo) and the result in fl2^ . 
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B.2 Robertson- Walker space-time 



We first discuss the case > 0. Similarly as in Appendix IB.ll it is convenient to use 
normalized vectors e^, = 0, l,r, that allow to write the embedding in (|29|1 as 



Z = f (a;) = 



1 



6(t)eo + a(t) C(r)ei + ^(r)e, 



for X = {t, r, 6,(p). A normal vector to all tangent vectors and such that h^^ = 1 is 

d{t) 



n4 



eo + 6(^) C(r)ei + 5(r)e, . 



By straightforward computations we find the following non-zero components of 



a 



40 



-J, 



a 



42 



a 



43 



b 

a 



and the result in (1531) . 

The computation for < is similar, and the formulas above remain true with the 
following changes: Bq and ei are exchanged, and VK is replaced by -\/|^- 

For K = we find 



lr/r2 



n4 = - 
2 



- + ro)d{t)-B{t) 
ro 



eo + - 



lr/r2 



ro )d{t) — B{t) ei + ra{t)er 

ro 



and a^^ as above with h = a. 



B.3 Early universe models 

We give some details on how we obtained the results reported in Section 15.21 We used 
MAPLE to perform the computations described below. 

We denote the l.h.s. of fl38|) and fl39|) as F and G, respectively. Inserting the ansatz in 
fj42|) we compute F = t^^^ (Fq + Fit + F2t^ + ■ ■ ■ ) for some sp (determined by x and y), 
and similarly for G. We then solve F„ = Gn = for n = 0, 1, 2, ■ ■ ■ . For low values of n 
(often = 0) we obtain conditions allowing us to determine x, y and the free parameters, 
and solving the equations for larger values of n gives equations allowing to compute 0„ 
and Hn in terms of the free parameters recursively. For the solutions thus obtained we 
compute the conservation law I in (140|) using fHTj) . This being constant is a useful check 
of our computations, and requiring this constant to be zero gives (usuall\l^ an additional 
constraint on the free parameters of our solutions. 

In the standard case x = we find three different solutions of (l38l) and (!39l) : (i) 
X = 4/3, ?/ = 1, 00 = a/2/A; no free parameter, (ii) x = 2/3, y = 0; one free parameter 

^■^Two exceptions are discussed in the last paragraph of this section. 
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00, (iii) X = y = 0; three free parameters (j)o, (pi, and Hq. Computing / we find that the 
solutions (i) and (ii) are incompatible with / = 0, and they therefore have to be discarded. 
For the solution (iii) the constraint / = determines Hq in terms of 0o and (pi, and we 
thus obtain the result given in f H5]) . For completeness we also give the coefficients of the 
0(t^)-terms of the solution (iii) in 



^,=l{m' + X4)Ho(l>o + ^ 



+ 



31 



[m" + \(t>l)HQ(, 



X 



AMI, 



+ 



PI 

001 
2 



A 



ml, 



^0 + 



m 



ml, 



A 
2J 



m 
~6 



2 A 

+ 77 



{m' + \(gf(P, 



ml, 



m 

h A 

6 



(59) 



In the extended case x = 1 we obtain six different solutions: (i) x = 3/4, y = 1, 
00 = l/(2v^2A); no free parameter, (n) x = 4/3, y = 1, (po = 1/(6\/A); no free parameter, 

(iii) X = (l + y3)/4 = 0.683(0), y = 0, 00 = ^6(^3 - 1)/(21 - liy3)Mpz = 1.50(2)Mpi; 
no free parameter, (iv) x = 1, y = 1/2; one free parameter 0o, (v) x = y = 0; five free 
parameters 0o, (pi, Hq, Hi, and H2, (vi) x = 1, y = 1/2; one free parameter 0o 7^ 0. 
Computing I we find that the solutions (ii), (iv) and (vi) are incompatible with 1 = 0, 
and they therefore have to be discarded. For the solutions (i) and (iii) we find that X 
in pT]) is identically zero, and these solutions are therefore allowed. By straightforward 
computations we find the results in flU]) and fllS]) for the solutions (i) and (iii), respectively. 
As discussed in Section IHT^ we believe that the solution (iii) in fH^ is of particular interest, 
and we therefore give the coefficients of the 0(t^)-terms in (145|) : 

03 =Mpi (^0.103(3) X'^Mj., + 0.000278(3) A^ + 0.00989(0) 

+ 0.00618(6) m^A + 0.0286(9) A ml, + 0.0766(3) m^XMl,^ (^gg^, 

ifg =0.00913(8) m^A + 0.00298(1) m^XMl, + 0.0109(9) X^M^, 
+ 0.00607(8) A^ + 0.0118(0) A XMl, - 0.000568(1) 

For the solution (v) we compute I in fHU]) using fHT]) . and we check that this is a constant 
depending on the free parameters of this solution. Demanding J = determines the free 
parameter H2 in terms of the other ones, and we thus obtain the solution in ( H6|) . It is 
straightforward to compute the coefficients of the 0(t'^)-terms in ( H6|) . but the result is 
lengthy, and we therefore do not give it here. 

The solution in (150|) -( |5T]) can be obtained by the same methods as the other solutions 
described above, using the ansatz in (H2|) with x = y = (pQ = 0. Alternatively, one can 
obtain this solution by a suitable limit from ( H6|) (the result for Hq is easy to understand 
from ( H^ . but the result for Hi is somewhat more subtle; note that the naive guess 
Hi = 2Hq is incorrect). 

It is interesting to note that, if one finds a solution of ( |38|) and ( 139|) with x 7^ and 2y 
an integer, this solution can have non-zero / in ( I40p only if 3x is an integer. Indeed, (141 p 
shows that X is a Laurent series in t, and, since a(tY is equal to t^^ times a Laurent series 
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in t (this follows from the definition of H{t) and fj42l) ). / is equal to t^^ times a Laurent 
series in t. Thus, if 3x is not an integer, / can be independent of t only if it is identically 
zero, and this is the case for our solutions (i) and (iii) in the case x = 1- For all other 
solutions of fl55]) and fl5^ with x ^ 0, 3x is an integer, and J is a constant depending on 
the free parameters. It is remarkable that, in all these latter cases, we find that J = is 
incompatible with other constraints. 



B.4 Higher-derivative Lagrangian systems: a toy model 

We consider a generalization of the harmonic oscillator defined by the following Lagrangian 

Lo = ^(q'-uV-af) (61) 

with a; > and a > constants and q = q{t) the dynamical variable. The Euler-Lagrange 
equation of this system is the following 4-th order differential equation, q + u'^q + a 4=0, 
and it has the following general solution 

q{t) = ci cos((X'+t) + C2 sin((X'+t) + C3 cos(w_t) + C4 sin(u;_t) (62) 
with four free parameters Cj and u± = ([1 ± (1 — 4:au'^y^'^]/{2a)Y^'^ , i.e., 

u_ = u + lu^a + 0{a'^), u+ = ^ - l^/^u^ + 0{a^^^). (63) 
2 y a 2 

We thus see that the limit a — )■ 0"*" of this system is delicate: the general solution describes 
a "slow" motion with angular frequency u^, which is a slightly perturbed motion of the 
harmonic oscillator (i.e., the solution of the model with a = 0, which only depends on 
two free parameters C1.2), but on top of it there is a "fast" oscillatory motion with angular 
frequency a;+ which diverges as a — ?► 0^. Thus the limit a — 0"^ in the model is only 
well-defined for solutions where this "fast" motion is absent, i.e., if 03 = 04 = 0. We 
can restrict to this subset of solutions which have a well-defined limit by the following 
constraints on initial conditions, 

q + colq = 0, q + ulq = at t = 0. (64) 
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